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The variational local moment approach (VLMA) solution of the single impurity Anderson model is 
presented. It generalizes the local moment approach of Logan et al. by invoking the variational principle 
to determine the lengths of local moments and orbital occupancies. We show that VLMA is a comprehen- 
sive, conserving and thermodynamically consistent approximation and treats both Fermi and non-Fermi 
liquid regimes as well as local moment phases on equal footing. We tested VLMA on selected problems. 
We solved the single- and multi-orbital impurity Anderson model in various regions of parameters, where 
different types of Kondo effects occur. The application of VLMA as an impurity solver of the dynamical 
mean-field theory, used to solve the multi-orbital Hubbard model, is also addressed. 
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I. INTRODUCTION 

The minimum at low temperatures in the resistance of 
metals with magnetic impurities^ the so-called Kondo 
effect, is related to generic singularities in the many- 
body perturbation expansion. 2 Many efforts have been 
put to developing a complete theory of the Kondo ef- 
fect, which would be free from such infinities 4^^^ 
Below a typical crossover temperature (Kondo temper- 
ature) there is a many-body screening of the localized 
magnetic moment by extended electrons with spini&i^ As 
a result, a many-body singlet state is formed. The local 
moment is shielded or confined and the properties of such 
a many-body system are well represented by the local 
Fermi liquid pictured This is a strong-coupling regime. 
Above the Kondo temperature the localized magnetic 
moment is decoupled from the electrons giving a Curie 
contribution to the magnetic susceptibility. This is a lo- 
cal moment regime. Now, the Kondo effect attracts at- 
tention elsewhere)^ in classica l 13 ! 14 and quantum phase 
transitions^ 5 in quantum dots and nanostructures^ in 
strongly correlated electrons, 17 and also in quark confine- 
ment and asymptotic freedom in hadronic matter j 18 ! 19 ! 20 

Traditionally, the Kondo effect is studied either within 
the Kondo model, 2 where the extended electrons and the 
localized magnetic moment interact via the exchange an- 
tiferromagnetic coupling, or within the single impurity 
Anderson model (SIAM), 3 where the extended electrons 
hybridize with the localized electrons and the latter in- 
teract via on-site Coulomb-like interaction. 

The Kondo model and the SIAM with linear dispersion 
relations have been solved exactly by the Bethe ansatz 
technique. 7,8 This powerful method is not, however, ap- 
plicable in general cases. Other analytical, asymptot- 
ically exact methods, i.e. conformal field theory and 
bosonization, are also limited to specific dispersion re- 
lations. In many practical cases, e.g. in modeling trans- 



port in nanostructuresi£ or in modeling correlated elec- 
tron systems within the dynamical mean-field theory 
(DMFT), 17 the Kondo model or the SIAM would have 
to be solved for arbitrary dispersion relations, exchange 
interaction, and/or hybridization functions. This re- 
quirement is achieved within numerically exact methods, 
e.g. in quantum Monte Carlo (QMC) 21,22 or in numeri- 
cal renormalization group (NRG)i&2£ The computational 
cost in time and memory increases, however, very fast at 
lowest temperatures or in systems with many channels 
and orbitals. Therefore, all those exact methods are of 
limited use in modeling real systems. 

Different analytical, approximate methods have also 
been developed. In fact, our understanding of 
the Kondo effect arrived from renormalization group 
approaches j 24 i 25 i 26 where high energy modes are succes- 
sively removed and low energy couplings are renormalized 
accordingly. Another class of approximate methods are 
those based on different perturbation expansions ! 27 ! 28 ! 29 
They do not reproduce either the exponential Kondo en- 
ergy scale or the local Fermi liquid properties at ener- 
gies below the Kondo energy scale. These drawbacks can 
be partially cured within sophisticated resummations of 
Feynman diagrams ! 30 ! 31 

In recent years a local moment approach 
(LMA) has been put forward by D. Logan and 
collaborators j 32 i 33 i 34 i 35 Within this analytical, approxi- 
mate method applied to the SIAM one recovers: i) the 
exponential Kondo energy scale, ii) low-energy local 
Fermi liquid properties, and iii) high energy charge 
fluctuations, correctly. The LMA starts with an unre- 
stricted Hartree-Fock (UHF) perturbation expansion, 
derived for each direction of the local moment, and 
afterwards the symmetry restoration is performed ] 32 ! 36 
Physically relevant observables are derived from aver- 
aging perturbation results over different local moment 
orientations. Finally, a free parameter, i.e. the length 
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of the local moment, is determined via imposing the 
local Fermi liquid conditions on the ground state. Since 
the LMA has been proved to reproduce several exactly 
known results and to recover correctly the properties of 
the SIAM at all energy scales^ 3 - 3 -* 3 ^ 3 ^ it seems to be 
an excellent method of choice in studying properties of 
nanosystems 16 and systems with correlated electrons^ 
Indeed the LMA has been used in studying Kondo in- 
sulators and heavy- fermions j 38 i 39 i 40 i 41 impurity systems 
with many orbitals, 42 and SIAM in the presence of 
magnetic field. 43 

In this paper we generalize the LMA of D. Logan et 
a/ j 32 i 33 i 34 i 35 by using the variational principle in deter- 
mining the length of the local magnetic moment. Such 
approach is motivated by the very general principle that 
the best approximate solution must minimize the ground 
state (free) energy. It is free from any bias toward the 
local Fermi liquid properties. Therefore, the variational 
local moment approach (VLMA) is applicable to any sys- 
tem with a magnetic impurity, not necessarily that be- 
ing in the local Fermi liquid fixed point. In particular, 
VLMA can be used in studying Kondo effects in multi- 
orbital and/or multi-channel SIAM, where local Fermi or 
non-Fermi liquid ground states are expected. 44 It can be 
extended to finite temperatures, where local Fermi liquid 
properties are only approximately satisfied. VLMA can 
also be applied to solving the impurity problems inside 
the DMFT and its cluster extensions i 45 i 46 i 47 This allows 
for unbiased studying of metal-insulator transitions in 
correlated electron systems in various configurations. In 
this paper we present a few applications of the VLMA. 

In the following we define the multi-orbital SIAM in 
Sec. II and next formulate the VLMA method in Sec. III. 
Here we also provide arguments that VLMA is a conserv- 
ing and thermodynamically consistent approximation in 
the Kadanoff-Baym sensed. To show the usefulness of 
this method, the VLMA is applied to one- and two- 
orbital SIAM in Sec. IV and to correlated electron prob- 
lems within the DMFT in Sec. V. Final conclusions and 
outlooks are in Sec. VI. 



II. MULTI-ORBITAL ANDERSON MODEL 



where S aG creates an impurity electron with spin a =|, \ 
and orbital quantum number a = 1,2, ...,iV, n aa = 
di, a d a(T is the number operator of impurity electrons, 
and c^ aa creates a conduction electron with spin a in a 
band a. We assume the same number of orbitals (bands) 
for local and conduction electrons. The Hamiltonian (pQ) 
contains local interactions U a between impurity electrons 
on the same orbital and local interactions U' a/3 between 
impurity electrons on different orbitals. The Ising part of 
the exchange interaction J a p is also included. The single 
particle energies of conduction electrons ek and the hy- 
bridization between impurity and conduction electrons, 
represented by matrix elements Vka/3 in (pQ) , can be com- 
bined into a hybridization function A(icj n ). It is a matrix 
in orbital space with matrix elements given by 

[A] a(3 (iu; n ) = A a/3 (i(j n ) = V . k , „ k , (2) 

^ lUJ n + fJL - 6 k 

where uo n are the fermionic Matsubara frequencies with 
(3 = l/(/c#T), and ksT is temperature in energy units, 
whereas jl is the chemical potential. 



III. VARIATIONAL LOCAL MOMENT 
APPROACH 

Starting from the LMA of Logan et al^^^^ we for- 
mulate its variational version, i.e. the VLMA, in order to 
approximately solve the Hamiltonian (pQ) with arbitrary 
parameters. The VLMA consists of four steps: 

i) finding the unrestricted Hartree-Fock (UHF) solu- 
tion of (pQ), 

ii) performing a renormalized perturbation expansion 
around each of the UHF solutions, 

iii) restoring the symmetry of the solution, 

iv) variationally determining the lengths of local mo- 
ments and orbital occupancies. 

Next we discuss each step in detail. 



The VLMA is developed explicitly for solving a multi- 
orbital single impurity Anderson model, which is defined 
first. The Hamiltonian of the multi-orbital SIAM is given 

by 

kaa kcra/3 




2^2 Yl ( U aP ~ Jap 5 **') navnpv', (1) 

<j<j' a^(3 



A. Unrestricted Hartree-Fock approximation 

The VLMA begins with the UHF solution of the SIAM 
in Eq. (pQ). To set the notation and explain basic ideas 
we first consider the one-orbital case, i.e. a = (3 = 1, 
U a = U, e a = e\. Then the UHF local Green's function 
G a (i(jj n ) (in the Matsubara representation) is given by 3 



lUJ n + /i - 61 - A(^ n ) - U (n-tr) 

where A(ico n ) = An(ico n ) in Eq. ([2]) and (n a ) is the 
average occupation number per spin. By introducing the 
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total average occupancy of the impurity n = (n^) + (rij) 
and the local magnetic moment 



/i = (n T ) - (n;), 
we write the UHF solution in the form 

G a (iuj n ) = - 



(3) 



with a = ±1 for spin up/down. The UHF solution de- 
pends on the value of the local moment \i = fiHF, which 
has to be determined self-consistently. It turns out that 
there exists a critical value U c such that for U > U c 
the UHF local moment \iuf ^ 0. 3 The sign (direction 
along the quantization axis) of the local moment is ei- 
ther positive or negative. The choice of one from the two 
possibilities breaks the symmetry of the solution with re- 
spect to the symmetry of the SIAM. This distinguishes 
two broken symmetry solutions, which we label by A for 
\i > and by B for \i < 0, i.e. 



G A J B ^n) 



1 



j n + ji - ei - A(icj n ) - If/n =F cr h u \f jL \ ' 



The Dyson equation 49 



(4) 



where G°(ico n ) 1 = iuo n + ft — e\ — A(ico n ) is the inverse 
of the non-interacting Green's function, defines two UHF 
self-energies Y^^ B = Un/2 ± aU/i/2. 

In the multi-orbital SIAM with N orbitals the UHF 
solution depends on 2N parameters: N local moments 
/jL a = (n a |) — (n a i) and N occupancies of the orbitals 
n a = (n a |)-h (n a |). In the UHF approximation the local 
moments have non-zero values above critical values of the 
interactions^ Similarly to the one-orbital case, the UHF 
solution with non-zero local moments breaks the symme- 
try of the Hamiltonian. In the case of non-degenerate or- 
bitals and zero magnetic field the UHF solution is doubly 
degenerate, i.e. solutions with {fi a } = {^l? M2, • • • , Mtv} 
and {— where all signs are opposite, have the same 
energy. These two equivalent solutions we denote again 
with A and £?, respectively. Hence, there are two Green's 
functions Ga' B (iu n ), which are now matrices in the or- 
bital space. The matrix Dyson equation 



t A/B 



defines the two matrix self-energies H^ B . The matrix 
elements of the self-energies are given by 50 

X)(2E£ 7 -J)r^^ (5) 

where \i = ^ a ja a is the length of the total impurity mag- 
netic moment. The UHF self-energy is diagonal in the 



orbital index a. Therefore, if the hybridization function 
([2]) is diagonal, the UHF propagators are also diagonal 
in the orbital indices, i.e. 

Only this diagonal case is considered here. 



B. Symmetry restoration 

Before discussing technical details of the perturbation 
expansion we first explain the idea of symmetry restora- 
tion. Since there are always two equivalent UHF so- 
lutions, the higher order in U a p and J a p perturbation 
expansion around the mean-field solution is not self- 
evident. We perform the renormalized perturbation ex- 
pansion (RPE) around each of the UHF solutions and 
then restore the symmetry. Following D. Logan ffi 33 i 34 i 35 
we formulate a symmetry restoring ansatz in the multi- 
orbital case as follows 

G a (i(j n ; {fi a }) = - (G^(zcj n ; {fi a }) + G B (iu n ; . 

Note, that the two RPE local Green's functions and 
G B are related by G^(iuj n ; {fi a }) = G B (iu n ; {- fi a }) . 
The result (j6]) depends on the values of all local moments 
H a and, in the general case, also on all local orbital occu- 
pancies n a . These are parameters of the method, which 
are not equal to their corresponding UHF values. These 
parameters are determined here via the variational prin- 
ciple (see Sec. HID). 



C. Renormalized perturbation expansion 

A IB 

Now, we need to find G a (ioj n ). We note that the 
symmetry restoration step does not depend on a spe- 
cific RPE scheme around UHF. However, it is important 
to select such an approximation that captures the pro- 
cesses which are essential to the Kondo problem. These 
are processes with the local spin-flips. 2,3 In the multi- 
orbital case spin-flip processes can occur together with 
the change of the orbital index. Such an orbital-flip is 
also possible without spin-flip, however. The second or- 
der Feynman diagrams describing these processes are de- 
picted in Fig. [TJ 

The second order contribution to the self-energy is 
given by 

(7) 

where v m is the bosonic Matsubara frequency and U the 
local interaction matrix with elements that have spin and 
orbital indices and are given by 



TTcrcr 

U a(3 



dapda-a'Ua + (1 - ^aa'J)- (8) 
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(a) c*Lo 






(c) c^c 




where we replaced °U AA (ih' n ) by n AA (i^ n ). 



FIG. 1: 2 nd order Feynman diagrams : (a) a spin- flip process 
without the change of orbital number; (b) spin- and orbital- 
flip; (c) only orbital- flip. The wavy lines correspond to the 
local interactions: (a) U a , (b) U^p, (c) U' a p — J a p- These 
diagrams contribute to the 2 nd order self-energy S^^(ia; n ). 



n= 



+ 



+ + 




(ico n ) have the sym- 



The UHF Green's functions G A/B 
metry 

G A (icu n ) = G* a (iu, n ). 

This leads to the relations between matrix elements of 
U BB (iu n ) and H. AA (iv n ), i.e. 

n^™\iv n ) = n A £-°-°\iv n ). (12) 

Hence, we have a useful relation between the two self- 
energies 

Y% a {iw n ) = Xi_ a (iuj n ). (13) 

The frequency-dependent self-energies S^( B (zo; n ) can 
now be used in the Dyson equation to obtain the two 

A IB 

propagators G a (i^n)- The full solution is given by 
the symmetry restoration equation ([6]). The solution is 
still dependent on the local moments and on the 

orbital occupancies. They are determined next. 

D. Variational principle 



FIG. 2: The ladder diagrams contributing to the full polar- 
ization propagator U AA (iv n )). 



By construction, since A(icu n ) is diagonal, the self-energy 
is also diagonal in orbital and spin indices. The polar- 
ization propagator °U^°' ° ' (iv m ) that contributes to the 
self-energy in Eq. ([7j) is given by 

°U A p ~ aa (ivn) = -77 ^ G A ~° '(i(jJm)G Ao \iv n + ioo m ), (9) 

iUJm 

in the transverse spin channel (a' = — cr), and by 



Ha/3 (IVn 



in the longitudinal channel (a' = a). Feynman diagrams 
representing these expressions are plotted in Fig[TJ 

Finite order perturbation theory does not reproduce 
the Kondo exponential scale. Therefore, following 
Refs. 32.33.34,35, we perform the ladder summation of 
diagrams with spin- and orbital-flip processes to infinite 
order. Then we obtain the self-energy with a new polar- 
ization propagator U AA (iv n ). The summed series, de- 
picted symbolically in Fig. [2j leads to a matrix expression 
for the polarization propagator 



n AA (iu n ) = °n^K)(i - \j°n AA {iv n ))-\ (10) 

where U is the interaction matrix (|5J). Then the self- 
energy is obtained by inserting the full polarization prop- 
agator U AA (iv n ) to the expression ©, i.e. 

T, A a (iuj n ) ^2(UZp ) 2 ^ ^2& Aa (iuj n -iiy m )Il AAc7 a (%Vm), 



The values of the local moments are yet unknown pa- 
rameters. The physical values are found by the minimiza- 
tion of the ground state energy at zero temperature or 
the free energy at finite temperatures. This variational 
approach is the main difference between the VLMA, pre- 
sented here, and the LMA in Refs. ^SUll • The 
general variational approach allows us to use the method 
for more than one orbital in the degenerate and non- 
degenerate cases, and also enables us to formulate the 
method directly at finite temperatures j 51 i 52 i 53 

In the SI AM, the ground state energy is obtained from 
the one-particle local Green's function. The total ground 
state energy in the SI AM has two contributions: i^buik 
and i^imp- The energy £?buik does not depend on the 
local moments {fi a }- Therefore, we only need to find the 
minimum of -Ei mp , which is given by 54 



2iri — j r 



G%a(«>), (14) 



2 



(11) 



where C is a contour in the complex frequency upper half- 
plane, and Gp a (u) and A a p(uj) are the retarded Green's 
and hybridization functions, respectively. After calculat- 
ing the integral, we express the ground state energy via 
real functions of real frequencies, i.e. 

Simp = (w + £a)0G2a(w) + 

1 ^2 F*± ^A a/3 (w)«G^ a (w)+«A a(S (w)9fG? a (w)) - 

Llj^ — a3 — ^h+u, — — — QGp^ui^iS) 
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The physical values of {fi a } and {n a } are found within 
the variational principle by minimizing Ei mp with respect 
to those parameters. 

It is important to observe, that the VLMA method can 
be implemented by using functions on real frequency axis. 
This gives us a direct access to the local spectral func- 
tions A^Juj) = — ^G^p(uj -f z0 + ) without the necessity 
of performing analytical continuation. The related cal- 
culations are presented in the Appendix. 

In order to use the VLMA at finite temperatures the 
free energy has to be minimized instead of the ground 
state energy. The calculation of the free energy from one- 
particle Green's function involves integration from zero 
to one over an auxiliary parameter A which multiplies 
the local interaction Since A enters in a nonlinear 
and indirect way into the expression for the free energy, 
the VLMA at finite temperatures is computationally in- 
volved and has not been implemented so far. The results 
presented in sections IV and V have been obtained at 
T = 0. 



E. Luttinger-Ward functional 

Any reliable approximation in the many-body theory 
should be thermodynamically consistent and conserving. 
I.e. physical observables should be the same irrespective 
of the way how they are determined and the microscopic 
conservation laws (e.g., conservation of energy, charge, 
spin, etc.) should be fulfilled. Within the Kadanoff and 
Baym formalism, 48 we formally show that the VLMA is 
a thermodynamically consistent and conserving approxi- 
mation. Its results are therefore reliable. In order to show 
this, it is enough to explicitly construct the Luttinger- 
Ward functional $[G] for VLMA. Here, $[G], which ful- 
fills the necessary condition that 6$[G]/5G = XI is 
given by 

$[G] = $[G A , G B ] = X - (4>[G A ] + 4>[G B ]) + 
±Trlog(G A G*) - Trlog Q (G A + G*)) , 

where Tr = T ^2 iUn tr, and the symbol tr means the trace 
in the orbital and spin space, and the functional <I> is the 
Luttinger-Ward functional in the ladder approximation 
used in Sec. IIIC. 55 It is constructed diagrammatically 
from G A and G B , respectively. The symmetry restoring 
constraint, G = |(G A + G S ), is also imposed. The 
free energy functional Q[G] is now expressed via 3>[G] 
and has the following form 

Q[G] = $[G] + Trlog(-G) - TrEG. 

The stationarity condition c5fi[G]/c5G = gives then the 
Dyson equation. 
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FIG. 3: (Color online). Local spectral functions for one- 
orbital SIAM with semi-elliptic hybridization (half- width D — 
10, n — 1) and different values of the interaction U. Inset: 
the energy regime around Fermi level at uj = 0. 



IV. RESULTS FOR ONE-ORBITAL SIAM 

Firstly, we test VLMA in the case of one-orbital SIAM. 
The spectral functions for the symmetric SIAM (n = 1) 
obtained using VLMA are presented in Fig. [3j Different 
spectral functions in this figure correspond to different 
values of the local interaction U a = U\ = U . The re- 
sults presented here are for a semi-elliptic hybridization 
function An (a;) = A (a;) defined by ([2]), whose imaginary 
part is given by 

^AH = (- A ^ M< c (16) 
1 M >D. 

where D is the so-called half-width and the real part 
5ft A (a; ) follows from the Kramers-Kronig relations. 



A. Symmetric one-orbital SIAM 

For U = the VLMA reproduces the exact result. The 
solution for U ^ has local Fermi liquid properties with 
the local spectral function pinned to its non-interacting 
value at the Fermi energy (cf. Fig. [3]) according to the 
Luttinger theorem^ As U increases, the characteris- 
tic three-peak structure emerges with the quasi-particle 
peak narrowing exponentially with U. The satellite peaks 
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^HF 




FIG. 4: (Color online). The value of the local moment as 
a function of the interaction U. The grey (red) curve corre- 
sponds to the VLMA solution (the minimum of Ei mp from 
Eq. (dH ). Black curve represents the self-consistent UHF so- 
lution \±hf- 



are the Hubbard bands. 57 Since the VLMA predicts the 
local Fermi liquid behavior, our results are in agreement 
with those obtained by the LMA if the LMA was applied 
to the hybridization function ([T6]) . 

In the symmetric SIAM with one orbital there is only 
one parameter, i.e. the value of the local moment /i, 
with respect to which the ground-state energy has to be 
minimized. The obtained value is plotted in Fig. [H as 
a function of the interaction U. For small interaction, 
the value of the local moment is equal to zero both in 
the VLMA solution (/iy in Fig. [J) and the UHF solution 
(fiHF in Fig. With the increase of U the values of 
local moments fly and [Ihf become different from zero 
at different J7, see Fig. 3J For large U the two solutions 
become exponentially close (with jiy always larger than 
Vhf)- 

We note that the closeness of the variationally obtained 
\iy to the UHF value \±hf makes searching for the mini- 
mum of the ground- state energy difficult at large U. This 
is due to the divergence of II (a;) in Eq. ([TO)) exactly when 
M = Vhf- Therefore, the numerical minimization of the 
ground state energy needs to be done carefully, especially 
for large U. 



B. Asymmetric one-orbital SIAM 

The results for asymmetric (n / 1) one-orbital SIAM 
are shown in Fig. [5j For fillings n < 1 the spectral weight 
is redistributed and, when n decreases, the Kondo reso- 
nance merges with the lower Hubbard band. The values 
of the spectral functions at the Fermi energy are still 
the same as those given by the Luttinger theorem^ The 
characteristic crossing points^ are also visible in Fig. [5l 

To conclude, the VLMA predicts the local Fermi liquid 
behavior as the solution with the lowest energy inside 
the manifold of different possibilities, which remains in 
accord to standard knowledge on the Kondo problem^ 
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FIG. 5: (Color online). Local spectral functions for different 
impurity occupancies. All plots for (7 = 4 and semi-elliptic 
hybridization with D — 10. The Fermi energy is at uj — 0. 



V. RESULTS FOR TWO-ORBITAL SIAM 

The two-orbital SIAM is physically richer than its one- 
orbital counterpart. Here the interplay between direct 
U a = U, U f a p = U' and exchange J a p = J interac- 
tions and the energy difference between various impurity 
orbital levels are of crucial importance. In certain pa- 
rameter regimes the Kondo effect is suppressed due to 
the quelling of fluctuations of spin and orbital degrees of 
freedom^ The VLMA confirms this prediction. 



A. SU(4) symmetric Kondo effect 

The results for the two-orbital degenerate case (ei = 
e 2 ) with SU(4) 44 symmetry (U = U' , J = 0) are pre- 
sented in Fig. El The local spectral functions have a three 
peak structure with the Kondo peak, which is pinned to 
the non-interacting value. The system remains in a local 
Fermi liquid fixed point up to arbitrary large values of 
U. The width of the Kondo resonance decreases expo- 
nentially with the value of U as it is in the one-orbital 
case. However, the prefactor in the exponent, which is 
proportional to the number of orbitals, makes the nar- 
rowing of the Kondo peak slower with increasing U as 
compared to the one-orbital case. To see this we empha- 
size different energy scales used on horizontal axis in the 
insets to Figs. [3] and El 



B. Influence of the exchange interaction 

In order to understand the influence of the non-zero ex- 
change interaction (Hund's coupling) we start from the 
atomic limit of the two-orbital SIAM (Vka/3 = 0)- I n this 
limit at J = and U = U' the impurity ground state 
with two electrons is six-fold degenerate. The presence 
of the Hund's coupling (J > 0) lifts this degeneracy, be- 
cause the states with higher total spin (5 = 1) have lower 
energy. The atomic ground state is in this case two-fold 
degenerate (S z = ±1). 
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FIG. 6: (Color online). Local spectral functions for one of the 
orbitals in the half-filled (n = 2) two-orbital SIAM for J = 0, 
U — U' and e\ — 62. The hybridization function is semi- 
elliptic with half-width D — 10, identical for both orbitals. 
Inset: the energy regime around the Fermi level at uj = 0. 




strong coupling 



FIG. 8: (Color online). Local spectral functions for one of 
the orbitals in the two-orbital SIAM for U — 2 and different 
values of the exchange interaction J (the plots correspond to 
regime (I) on the diagram in Fig. 0. Other parameters as in 
Fig. [6] Inset: the energy regime around the Fermi level at 

LU = 0. 




FIG. 7: (Color online). Phase space diagram (U, J) of the 
ground state of the two-orbital particle- hole symmetric SIAM. 
The dotted (brown) lines are plotted to aid in finding points 
in this diagram that correspond to spectral functions plotted 
in Figs. [8] and [9] The blank area corresponds to U' < (see 
text). 



FIG. 9: (Color online). Local spectral functions for one of 
the orbitals in the two-orbital SIAM for U = 10 and different 
values of the exchange interaction J (the plots correspond 
to regimes (II) and (III) on the diagram in Fig. 0. Other 
parameters as in Fig. [6] Inset: the area around the Fermi 
energy at co = 0. 



With non-zero hybridization, the impurity spin can be 
flipped due to higher order processes in Vka/3- This leads 
to full or partial screening of the impurity magnetic mo- 
ment. The VLMA solution of SIAM enables us to dis- 
tinguish three parameter regimes with different ground 
states: weak coupling (I), strong coupling (II), and local 
moment (III) regimes. These regimes are presented on 
the ground state phase diagram in Fig. We consider 
here the symmetric case with U' = U — 2 J. Therefore, 
for J < 0.5/7, U' becomes negative and this case is not 
analyzed here (blank area in Fig. [7j). The spectral func- 
tions for selected parameters in the regimes (I)- (III) are 
presented in Figs. [8] and [9l 

In the weak coupling regime (area (I) in Fig. [7j) the 
spins of the localized electrons are screened individu- 
ally and the exchange interaction J weakly influences the 
width of the central peak (cf. Fig. [8]). The weak coupling 



regime coincides with the regime where the UHF solution 
gives = fi2 F = (note however, that the values of 
the local moments in the VLMA are finite here). When 
U increases and when the UHF solution gives non-zero 
values of fi^ F and the Kondo resonance narrows 

and its width strongly depends on J (cf. Fig. [9]). This 
is the strong coupling regime (II) in Fig. [3 In the local 
moment regime, (area (III) in Fig. [7j) the Kondo reso- 
nance is absent. Now, the electrons, which have only 
spin 1/2, cannot fully screen the spin 1 on the impurity. 
The disappearing of the Kondo resonance gives rise to a 
critical line J C (U) that separates the strong coupling (II) 
and local moment (III) regimes (cf. Fig. [7j)- 

Finally we note that the local moment regime is ab- 
sent in the two-band two-orbital SIAM with the Heisen- 
berg exchange interaction. 59 Hence, the validity of the 
diagram presented here is restricted to the Ising-type ex- 
change interaction. 
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FIG. 10: (Color online) Local spectral functions in case of 
degenerate Se = and non-degenerate Se = 2 two orbitals in 
the half- filled case n — 2 and [7 = 6, J — 0; Fermi energy is 
at uj — 0. 



C. Non-degenerate orbitals 

The case of non-degenerate orbitals (Se = ci — €2 7^ 
0) is numerically more involved because of the increased 
number of parameters with respect to which the ground 
state energy (fT4|) has to be minimized. These parameters 
are: two values of the orbital local moments ii a (a = 
1,2), and two occupancies n a of the orbitals. The total 
occupancy n of the impurity can be adjusted by changing 
the chemical potential. Here we fix n = 2 (half-filling in 
the two-orbital case). 

Once the degeneracy is removed (Se > J) in the half- 
filled case both impurity electrons occupy the lower level 
and the total spin of the impurity is zero. The resulting 
local spectral densities show no quasi-particle peak at the 
Fermi energy (cf. Fig. [T0|). Fluctuations of the spin and 
orbital degrees of freedom are suppressed and the Kondo 
effect is absent. 

In conclusion, the VLMA correctly provides approxi- 
mate solutions to one- and two-orbital SIAM. Our results 
are in agreement with others 9,44 which shows that the 
VLMA is a good method for solving analogous magnetic 
impurity problems. 

VI. MULTI-ORBITAL HUBBARD MODEL 

In this section we apply VLMA to solve the DMFT 
equation s 17 ! 37 ! 60 for multi-orbital systems. The Hamilto- 
nian of the multi-orbital Hubbard model is given by 61 

^Hubbard = Hkin + Hy = 

^ ^2 ^2 ( U '<*P ~ j <*P s <t<t') n ia(T n i(3(J f , (17) 

i,cr,(j' a^(3 

where the interaction part Hjj is local and identical to 
the interaction part of the Anderson impurity model. In 
the DMFT, the lattice model (in this case the Hubbard 
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FIG. 11: (Color online). Local spectral functions for the one- 
orbital Hubbard model at half-filling. Left figure: the re- 
distribution of the spectral weight with increasing U due to 
Mott-Hubbard metal-insulator transition. Right figure illus- 
trates hysteresis - the solutions are obtained from initially 
insulating density of states 

model) is mapped onto a multi-orbital impurity model 
with the same local interaction. The impurity is coupled 
to a bath of particles, whose properties are determined 
self-consistently. In order to solve the resulting equa- 
tions, the single impurity problem has to be solved iter- 
atively, in each step with a new hybridization function 
until self-consistency is reached. The VLMA can be used 
for arbitrary hybridization function and therefore it can 
serve as an impurity solver for the DMFT. 

A. Results for one and two-orbital Hubbard model 

The results of the DMFT+VLMA approach for the 
one-orbital Hubbard model on the Bethe lattice 6 -^ with 
half-bandwidth D = 1 are shown in Fig. [Til The pa- 
rameters here are: interaction U a =i = U, filling n = 1, 
Fermi energy at jl = and zero temperature. The 
VLMA solver for DMFT predicts a Mott-Hubbard metal- 
insulator transition 63 with U& 7^ U C 2 (hysteresis). In the 
current implementation of the method the value of U c \ 
is obtained rather precisely (U c \ = 1.45). In the case 
of ?7 C 2, due to the exponential narrowing of the quasi- 
particle peak, the minimum of the ground state energy is 
increasingly difficult to find. That makes U C 2 harder to 
determine. The obtained value U C 2 ^2.5 may be under- 
estimated due to the difficulty in finding the minimum 
\iy close to a divergence of the polarization propagator 
(as discussed in Sec. IV). 

We note that in the solution of the DMFT equations 
within the LMA the hysteresis was almost absent. 64 This 
is due to the imposed Fermi liquid condition inside the 
LMA, which favors the metallic state. The VLMA is not 
biased in this respect. 

The Mott-Hubbard metal-insulator transition in two- 
orbital case also shows hysteresis. In this case the tran- 
sition is strongly influenced by the presence of non-zero 
exchange interaction J a p = J. The exchange interaction 
suppresses the Kondo effect and causes the transition to 



9 




FIG. 12: (Color online). Local spectral functions for the two- 
orbital Hubbard model with degenerate orbitals. Left: the 
redistribution of the spectral weight with increasing U due to 
Mott-Hubbard metal-insulator transition with no Hund's rule 
coupling J = 0. Right: the metal-insulator transition in the 
case of J 7^ 



merically in this regime. The potential difficulty at finite 
temperatures lies in an efficient calculation of the free en- 
ergy, as mentioned in Sec. HID. The finite temperature 
numerical implementation may, however, be free from the 
numerical problem of finding the minimum close to a di- 
vergence, which we faced in the ground state case. These 
issues are also left for the nearest future. 
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take place at a smaller value of the interaction U a = U 
as shown in Fig. [12] (here: U' a p = U — 2 J). 

In the two-orbital case, the DMFT solutions of the 
Hubbard model within VLMA show a metal-insulator 
transition that is orbital selective ] 51 ! 52 Such a transition 
is present both in the different bandwidths 52 case and in 
the non-degenerate orbitals case^i 



VII. SUMMARY 

We have formulated in Sec. Ill the variational local mo- 
ment approach and have applied it in Sec. IV to solve dif- 
ferent versions of the single impurity Anderson problem. 
The VLMA is a comprehensive, thermodynamically con- 
sistent and conserving approximation, as we have shown 
in Sec. HIE. It can be used for any values of the input 
parameters. Such a comprehensive method for solving 
the SIAM is a good method of choice for an efficient 
impurity solver in DMFT. We have presented in Sec. V 
some results for the one- and two-orbital Hubbard model 
obtained within the DMFT with VLMA. 

The main point where the VLMA is different from the 
LMA of Logan et al2&2&2&& is that the lengths of local 
moments and orbital occupancies are determined varia- 
tionally. Our method treats on equal footing both Fermi 
and non-Fermi liquid ground states as well as insulating, 
local moment states. 

The subtlety, however, in finding the ground state 
energy minimum at large interaction values, in partic- 
ular close to the Mott-Hubbard metal-insulator transi- 
tion, implies that at the moment the VLMA cannot yet 
be used as a "black-box" impurity solver. The semi- 
analytical formulation and implementation on the real- 
frequency axis make the method, in our opinion, valuable 
and we hope that the numerical difficulties of the current 
implementation will be overcome in the nearest future. 

The method has been formulated for arbitrary tem- 
perature. However, it has not yet been implemented nu- 



APPENDIX: REAL FREQUENCY 
IMPLEMENTATION OF THE RPE USED IN 
SEC. IIIC 

We introduce the spectral representation Da^ Ba (ou) 
for the UHF Green's functions (|4|) in the following way 

/+oo A/ B a / i \ 
,U a [00 ) 
duo — 
-oo lUJ n ~ U 1 

where D^ Ba {uj) is thus given by 

D* /BcT (oo) = -l%G£ /BcT (oo + iO + ). 

Then the polarization propag ator °U^- aa (iu n ) from 

([9]) is expressed by the spectral functions D^ a (uu) as fol- 
lows 



O-qAA — <J<7 f A 



up 



(iVn)-- 



+OG duj'duj"D^{uj')D A -^'^ 



(iVn + iUrn ~ 00')(iuO m - UO") ' 



Performing the sum over Matsubara frequencies 49 we ob- 
tain 



O-i-rAA — aa / ■ 

U a(3 (iv, 



»-//. 



du duo — : X 



IV n — 00" 

(f(uo f )-f(uo f -uo")), 

with f{uj) = (1 + exp(/3c^)) _1 being the Fermi-Dirac dis- 
tribution function. If we introduce the spectral represen- 

AA 



tation °Xad acr i 00 ) a ^ so f° r ^ ne polarization propagator 



Ot-t AAcrcr' / • \ f 
H a/3 \lVn) I 



+ oo AA aa' ( I 



duo 



AAaa ( f\ 
I Xcx(3 \U_ ) 

iVn — UO 1 



with 



AAaa 1 / \ 1 c\0tt AA acT ' ( , - n 4 



we can express °Xap acr ( u ) D Y integrals over real fre- 
quencies with the UHF spectral functions in the following 

way 

/+oo 
-oo 
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We can analogously express the other (orbital) polar- 
ization propagator °II^ crcr (^ n ) by integrals over the 
UHF spectral functions. The real parts of the polariza- 
tion propagators can be found by using Kramers-Kronig 
relations^ 

The sum over the Matsubara frequencies in Eq. ([TT]) is 
done with the use of the spectral representations of the 
full polarization propagator Yl AA (ii/ n ), i.e. 



nAA era 1 / • \ f 
a/3 0*M = / 



+ 00 AAaa' ( I 



doo 



AAaa ( l\ 
/Xa/3 \ U ) 



00 ^ri 



where 



AA aa' / \ 1 r^-r-rAA aa' / , • rv 4 

X a( 3 [uj) = -- ^sU a(3 (oo + 

TT 



Then Eq. ([TT]) can be rewritten to give 

^ doo'doo"D^'{oo')x^ a ' a {u") 

v^a/3 ; 1 J 

{3 a' %Vr, 



(lUJn — iVm — Oj')(iVm ~ UO") 



After performing the sum over Matsubara frequencies we 
arrive at the expression 



RrT , jj-00 



„ pr (u" - w') 

duo — — x 



lUO n — UO 



which contains the Fermi-Dirac distribution f(uo) to- 
gether with a function f(oo) = (1 — e^) -1 . Introducing 
the spectral representation also for the self-energy 



^Aa f . \ [ +OC j B^ a (i0) 

L a iiuon) — I duo- 



with 



B^(uo) 



we are able to calculate the self-energy using only inte- 
grals over real frequencies 



. r+oc 

B^(oo) = ^2(U% f / dJDf {oo - oo') x 

(3a> J -°° 



Xf3a a a (oo')(f(oo' -oo)- f(oo')). 



The real part of self-energy can be obtained from the 
imaginary part with the use of Kramers-Kronig relations. 



* Electronic address: Anna.Kauch@physik.uni-augsburg.de 

1 W. Meissner and B. Voigt, Ann. Phys. 7, 761; ibid. 892 
(1930). 

2 J. Kondo, Prog. Theor. Phys. 32, 37 (1964). 

3 P. W. Anderson, Phys. Rev. 124, 11 (1961). 

4 J. R. Schrieffer and P. A. Wolff, Phys. Rev. 149, 491 
(1966). 

5 K. G. Wilson, Rev. Mod. Phys. 47, 773 (1975). 

6 H. R. Krishna-murthy, J. W. Wilkins, and K. G. Wilson, 
Phys. Rev. B 21, 1003 (1980). 

7 A. M. Tsvelik and P. Wiegmann, Adv. in Phys. 32, 453 
(1983). 

8 N. Andrei, K. Furuya, and J. H. Lowenstein, Rev. Mod. 
Phys. 55, 331 (1983). 

9 A. Hewson, The Kondo problem to heavy fermions (Cam- 
bridge University Press, 1993). 

10 P. Coleman, in Lectures on the physics of highly correlated 
electron systems, edited by F. Mancini (American Institute 
of Physics, 2002), vol. 629, pp. 79-160. 

11 P. Nozieres, Theory of interacting Fermi systems (West- 
view Press, 1997). 

12 G. A. Gehring, J. Phys. Condens Matter 14, V5 (2002). 

13 J. M. Kosterlitz and D. J. Thouless, Prog. Low Temp. 
Phys. VIIB, 371 (1978). 

14 A. O. Gogolin, A. A. Nersesyan, and A. M. Tsvelik, 
Bosonization and strongly correlated systems (Cambridge 
University Press, 1998). 

15 S. Sachdev, Quantum phase transitions (Cambridge Uni- 
versity Press, 2001). 

16 L. Kouwenhoven and L. Glazman, Physics World 14, 33 
(2001). 



17 G. Kotliar and D. Vollhardt, Physics Today 3, 53 (2004). 

18 K. G. Wilson, Rev. Mod. Phys. 55, 583 (1983). 

19 CM. Varma, in Windsurfing the Fermi Sea (Elsevier Sci- 
ence Publishers B.V., 1987), vol. 1, p. 69. 

20 P. Coleman, Ann. Henri Poincare 4, S559 (2003). 

21 R. Fye and J. Hirsch, Phys. Rev. B 38, 433 (1988). 

22 P. Werner, A. Comanac, L. de' Medici, M. Troyer, and 
A. J. Millis, Phys. Rev. Lett. 97, 076405 (2006). 

23 R. Bulla, T. A. Costi, and T. Pruschke, Rev. Mod. Phys. 
80, 395 (2008). 

24 P. W. Anderson, J. Phys. C: Solid State Phys. 3, 2436 
(1970). 

25 S. Kehrein, The flow equation approach to many-particle 
systems (Springer, 2006). 

26 C. Karrasch, R. Hedden, R. Peters, T. Pruschke, 
K. Schonhammer, and V. Meden, J. Phys.: Condensed 
Matter 20, 345205 (2008). 

27 V. Zlatic and B. Horvatic, Phys. Rev. B 28, 6904 (1983). 

28 K. Yosida and K. Yamada, Prog. Theor. Phys. Suppl. 46, 
244 (1970). 

29 T. A. Costi, J. Kroha, and P. Woffle, Phys. Rev. B 53, 
1850 (1996). 

30 V. Janis and P. Augustinsky, Phys. Rev. B 75, 165108 
(2007). 

31 I. S. Krivenko, A. N. Rubtsov, M. I. Katsnelson, and A. I. 
Lichtenstein, cond-mat/0910.0792. 

32 D. E. Logan, M. P. Eastwood, and M. A. Tusch, J. Phys. 
Condens. Matter 10, 2673 (1998). 

33 D. E. Logan and M. T. Glossop, J. Phys. Condens. Matter 
12, 985 (2000). 

34 D. E. Logan and N. L. Dickens, J. Phys. Condens. Matter 



11 



14, 3605 (2002). 

35 M. T. Glossop and D. E. Logan, J. Phys. Condens. Matter 
14, 6737 (2002). 

36 W. S. Verwoerd, Phys. Lett. A 43, 535 (1973). 

37 A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg, 
Rev. Mod. Phys. 68, 13 (1996). 

38 V. E. Smith, D. E. Logan, and H. R. Krishnamurthy, Eur. 
Phys. J. B 32, 49 (2003). 

39 N. S. Vidhyadhiraja, V. E. Smith, D. E. Logan, and H. R. 
Krishnamurthy, J. Phys. Condens. Matter 15, 4045 (2003). 

40 N. S. Vidhyadhiraja and D. E. Logan, Eur. Phys. J. B 39, 
313 (2004). 

41 D. E. Logan and N. S. Vidhyadhiraja, J. Phys. Condens. 
Matter 17, 2935 (2005). 

42 M. R. Galpin, A. B. Gilbert, and D. E. Logan, J. Phys. 
Condens. Matter 21, 375602 (2009). 

43 C. H. Kim and J. Yu, J. Phys. Condens. Matter 19, 456203 
(2007). 

44 D. L. Cox and A. Zawadowski, Adv. Phys. 47, 599 (1998). 

45 A. N. Rubtsov, M. I. Katsnelson, and A. I. Lichtenstein, 
Phys. Rev. B 77, 033101 (2008). 

46 T. Maier, M. Jarrell, T. Pruschke, and M. H. Hettler, Rev. 
Mod. Phys. 77, 1027 (2005). 

47 G. Kotliar, S. Y. Savrasov, G. Palsson, and G. Biroli, Phys. 
Rev. Lett. 87, 186401 (2001). 

48 G. Baym and L. KadanofT, Phys. Rev. 124, 287 (1961). 



49 A. I. Fetter and J. D. Walecka, Quantum theory of many- 
particle systems (McGraw Hill, 1988). 

50 P. B. Coqblin and A. Blandin, Adv. Phys. 17, 281 (1968). 

51 A. Kauch and K. Byczuk, Physica B 378-380, 297 (2006). 

52 A. Kauch and K. Byczuk, in Quantum magnetism, edited 
by B. Barbara, Y. Imry, G. Sawatzky, and P. C. E. Stamp 
(Springer, 2008), p. 79. 

53 A. Kauch, Ph.D. thesis, University of Warsaw (2009). 

54 B. Kjollerstrom, D. Scalapino, and J. SchriefTer, Phys. Rev. 
148, 665 (1966). 

55 J. W. Orland and H. Negele, Quantum many particle sys- 
tems (Westview Press, 1998). 

56 J. M. Luttinger, Phys. Rev. 119, 1153 (1960). 

57 T. Pruschke and N. Grewe, Z. Phys. 74, 439 (1989). 

58 M. Eckstein, M. Kollar, and D. Vollhardt, J. Low Temp. 
Phys. 147 (2007). 

59 T. Pruschke and R. Bulla, Eur. Phys. J. B 44, 217 (2005). 

60 W. Metzner and D. Vollhardt, Phys. Rev. Lett. 62, 324 
(1989). 

61 J. Hubbard, Proc. R. Soc. A281, 401 (1964). 

62 M. Eckstein, M. Kollar, K. Byczuk, and D. Vollhardt, 
Phys. Rev. B 71, 235119 (2005). 

63 N. F. Mott, Rev. Mod. Phys. 40, 677 (1968). 

64 M. P. Eastwood, Ph.D. thesis, Oxford University (1998). 



